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Abstract 

Correlation between shower partons is first studied in high px jets. Then in the 
framework of parton recombination the correlation between pions in heavy-ion collisions 
is investigated. Since thermal partons play very different roles in central and peripheral 
collisions, it is found that the correlation functions of the produced hadrons behave 
very differently at different centralities, especially at intermediate px- The correlation 
function that can best exhibit the distinctive features is suggested. There is not a 
great deal of overlap between what we can calculate and what has been measured. 
Nevertheless, some aspects of our results compare favorably with experimental data. 
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1 Introduction 


Correlations between hadrons produced in jets in heavy-ion collisions have become a sub¬ 
ject of intense study recently [1]-[5]. There are many ways to describe the nature of the 
correlations, some using triggers, others giving correlations without triggers. So far no clear 
understanding of the properties of dihadron correlations has emerged. Satisfactory theo¬ 
retical description of the problem is lagging even farther behind. In this paper we offer a 
theoretical view of some aspect of the correlations among partons and hadrons in connection 
with jets. Because we treat only a limited aspect of the dihadron distributions (notably only 
7r + 7r + correlation), we cannot make exhaustive comparisons with the data that have their 
own limitations, such as the lack of particle identification in many cases. Nevertheless, our 
results illuminate some important features of the nature of the correlations. 

For experiments at the Relativistic Heavy Ion Collider (RHIC) the relevant range of trans¬ 
verse momenta ( p ?) of the detected hadrons in which the correlations are probed is between 
0.5 and 9 GeV/c. Since no calculation from first principles is reliable in the intermediate pt 
range, any theoretical treatment of correlations can only be done in the framework of some 
phenomenological model. The parton recombination model has been successful in describing 
the single-particle inclusive distributions in Au+Au [6] as well as d+Au collisions [7, 8]. 
Indeed, it has already been applied to the study of some aspects of dihadron correlations [9]. 
We shall extend that study and examine in more detail various features of the correlations, 
not only among the produced hadrons, but also among the partons before they hadronize. 

There are limitations to what we can do here. We shall continue to use the one¬ 
dimensional (ID) formulation of parton recombination, which offers a transparent display in 
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analytical forms the various components that contribute to dihadron correlations, without 
any thing hidden in an elaborate code. Focusing only on the px variable in the ID treatment 
means that we cannot consider the dependencies on the pseudorapidity ( 77 ) and azimuthal 
angle (0) separations between the two detected hadrons. That shortcoming will be overcome 
in a separate investigation. But by not being distracted by the behavior in A?/ and A0 vari¬ 
ables, we can examine more closely the separate roles that thermal and shower partons play 
in the correlation problem, when only the collinear momenta are considered. Experimen¬ 
tally, they correspond to the magnitudes of the associated particle momenta. Even with that 
simplification we have not considered all charged particles, which are what the experiments 
measure without particle identification. To gain an overview of the general properties of the 
correlation functions, the hadron pair 7 r + 7 r + is sufficient and more transparent in providing 
a rich display of various features. 

Our basic assumption in this work is that the shower partons are emitted independently 
in a jet except for the momentum constraint that the sum of the momentum fractions of 
all partons cannot exceed 1. This assumption is only a working hypothesis in this study 
and is a reasonable starting point of our investigation of correlations in order to form a 
base line, from which the effect of dynamical correlation can be differentiated in the future. 
Whatever hadronic correlation that we shall calculate in this paper therefore follows from 
the kinematical correlation of the shower partons in a jet. The dynamical independence 
of the shower partons may perhaps be approximately correct at small momentum fractions 
where many shower partons are likely to be present. At larger momentum fractions the 
kinematical constraint is important and is suitably taken into account. Since the distributions 
of the shower partons have been determined from the fragmentation functions (FFs) with the 
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same assumption of dynamical independence [10], we are in that respect merely being self- 
consistent in our study of their correlation here. Since the single-particle distributions that 
we have calculated using those distributions of the shower partons are in good agreement 
with the data [6], there is indirect evidence that our assumption may not be too far from 
reality. 

Furthermore, another test of the reliability of the shower parton distributions (SPDs) 
has been carried out in [11], where the fragmentations of u and d quarks to p + p have been 
calculated using the SPDs for recombination into a proton or antiproton. The resultant FFs 
agree with the parametrization of KKP [12] over 4 orders of magnitude within a factor of 2, 
which is roughly the margin of error made in the KKP parametrization of the data on p + p 
production [11], Because of that error in [12], the work on the proton FF is not included 
in [10], but it is nevertheless a demonstration that the assumption of independent shower 
partons except for kinematical constraint can lead to sensible result in the calculation of 
the proton FF, which has not been done in any other approach to hadronization, as far 
as we know. In this paper we continue to use that assumption to calculate the hadronic 
correlations in heavy-ion collisions. 

We shall first examine the single-particle distributions for An+An collisions at all cen¬ 
tralities. The parameters determined can then be used for subsequent studies of correlations 
among partons and hadrons. In Sec. 3 the correlation among the shower partons is calculated 
both for jets in vacuum and for jets produced in heavy-ion collisions (HIC). They behave 
quite differently. Then we study in Sec. 4 the correlation among hadrons, for which the role 
of thermal partons becomes important in central collisions, but not as much in peripheral 
collisions. In Sec. 5 our results are compared with existing experimental data in areas where 


4 



such comparisons are possible. 


2 Single-Particle Distributions at all Centralities 

In [6] the single-particle distribution of hadrons produced in Au+Au collisions at 200 GeV 
is studied for 0-10% centrality only. We now complete the picture by considering all other 
centralities, thereby determining the corresponding parameters for the thermal partons and 
the value of £ that describes the average number of hard partons emerging from the dense 
medium to hadronize. This also gives us the opportunity to reiterate the basic equations for 
parton recombination. 

The invariant inclusive distribution for a produced meson with momentum p is [6] 


p ~T L = I ~~ F jj'(Qu < b) R *( < luq2,p) (1) 

dp J qi q 2 

where the recombination function (RF) for a pion is [13] 

= + ( 2 ) 

p \p p J 

The distribution for two partons, j and j', that are to recombine has three components 


Fjy = TT + TS + SS, (3) 

where the thermal parton distribution is 

T(q) = Cqe-' T , (4) 

and the shower parton in TS has the distribution 

<%) dkkfi(k)Sj(q/k) , (5) 
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whereas the two shower partons in the same jet (55) have the correlated distribution 


(S5)(«i, ®) = { £ / dkkf,(k) S{ 


<?2 


siM, si 


' ( x 2 


1 — X\ 


s 3 t(xi)si 


' / X 2 


1 — Xi 


+ 5/ 


k-qi / 
x 1 


1 - X 2 


Si (; X 2 ) 


( 6 ) 

( 7 ) 


The distribution _/)(&) of hard parton i at midrapidity in Au+Au collisions is given in [14], 
The integral over k is from k m i n = 3 GeV/c to an upper limit which we shall set at 40 GeV/c. 
Sj is a matrix of shower parton distributions (SPD) that have been determined in [10]. The 
form of the two SPDs, given in Eq. (7), symmetrizes the order of emission and guarantees 
the momentum conservation condition x\ + x 2 < 1, apart from which the shower partons 
are assumed to be independent. Thus the correlation between the partons is essentially 
kinematical, as it is treated here, although other correlations of more dynamical origin can 
be considered later when the need arises. 

The shower partons can either recombine among themselves, or recombine with the ther¬ 
mal partons. In the former case the process reproduces the FF 

xDJ(x) = f —— (si(x i), S{ R„(x 1 ,x 2 ,x) . (8) 

J x\ x 2 l V1 — X\) ) 


Indeed, this is the equation from which the SPDs are determined in the Erst place [10]. Thus 
the 55 term in Eq. (3), when substituted into Eq. (1), yields the fragmentation component 


P- 


d N* s 
dp 


= f£ I dkkf,(k){D’(‘j ), 


k 


V 


(9) 


Note that because the two shower partons are from the same jet, Eq. (9) depends linearly 
on ^/j(/c), which is the same dependence that the TS component has. Hence, when all three 
terms in and the RF R n are substituted into Eq. (1), we have, in the notation p = p T , 
dN« C 2 _„ /T <£ 


pdp 


= if^ /T + 1 E / l D ’ (f) + y / d i' s < (|) <p - 

IJ L 


. ( 10 ) 
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This is the basic formula for the inclusive distribution of pion. The hard parton i is summed 
over all relevant species, while j is summed over the two possible constituent quark species 
of the pion. The quark j picks up a partner from the thermal environment to form the pion. 
That TS component turns out to be dominant in the intermediate px region (3 < px < 9 
GeV/c) in central Au+Au collisions, and is responsible for the high p/n ratio, when the 
production of proton is considered in a similar way [6]. 

We now apply Eq. (10) to non-central Au+Au collisions. The hard parton distribution 
fi(k) is scaled down from central collisions by the average number of binary collisions (N co \\) • 
Our procedure is to fit the low-pr data (for p T < 2 GeV/c) by the first term that arises from 
TT recombination, thereby determining C and T for each centrality. The only remaining 
parameter, £, is then determined by the overall normalization of the spectrum at each cen¬ 
trality so that the shape of the large px dependence is a prediction of Eq. (10). The data 
used are the 7r° spectra from PHENIX that extend up to px = 10 GeV/c [15]. We, however, 
calculate the 7r + distribution with j = u and d, since our formalism does not distinguish 
the two. Figure 1 shows the result of our calculation for four representative centralities, 
although we have treated all nine centralities, for which the parameters are shown in Fig. 
2. Evidently, the spectra are well fitted at all centralities. The inverse slope T shows no 
significant dependence on centrality; the common value may be set at 

T = 0.305 GeV/c. (11) 

The value of £ increases from 0.09 at 0-10% to 0.97 at 80-92% centrality. Thus Au+Au 
collision at 80-92% is essentially equivalent to pp collision, since the suppression factor due to 
energy loss is nearly 1, i.e., no jet quenching. The value of C decreases from 21.8 (GeV/c)^ 1 
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at 0-10% to 2.3 (GeV/c) -1 at 80-92%. The thermal component in peripheral collision is 
therefore much reduced, as expected. 

The contributions from the various components are shown in Fig. 3 for 0-10% and 80-92% 
centralities. It is clear that at 0-10% the component from thermal-shower recombination is 
dominant in the 3 < pr < 9 GeV/c range, but at 80-92% that component is lower than 
both thermal-thermal and shower-shower components. Indeed, SS dominates for all pr > 3 
GeV/c; it means that the fragmentation model is quite adequate for pp collisions in that pt 
range. 


3 Correlation of Partons in Jets 


Before we consider the correlation between hadrons produced in jets, let us first study the 
simpler problem of correlation between partons in jets. Although the latter is not measurable, 
it provides a striking view of the difference between jets produced in vacuum, as in e + e“ 
annihilation, and jets produced in HIC where the hard parton momentum is not fixed. 

To define our notation, let pi(l) and p 2 ( 1,2) be the single-particle and two-particle dis¬ 
tributions, and we define the two-particle correlation function as 

C 2 (l, 2) = p 2 (l, 2) — p 1 (l)p 1 (2) . (12) 


Later, when we apply this formula to hadron correlation, the experimentally measurable 
quantities are 


Pi(l) 


dN ni 

Pidpi ’ 


P2( 1,2) 


dN n 


Pidpip 2 dp 2 ' 


(13) 



The normalized correlation function is 


K ^=Mtkr r ^ 2) - 1 


(14) 


where 


r 2 ( 1 , 2 ) = 


P2{ 1 , 2 ) 


pi(l)pi(2) 

For reasons that will become evident below, it is useful to consider the ratio 


(15) 


Gzi 1 ) 2 ) — 


C 2 (l,2) 


(16) 


[Pi(l)l>i(2)]‘ /2 ’ 

a quantity that is similar to the Pearson’s correlation [16] under certain assumptions, and is 
closely related to a quantity considered in [2] 


3.1 Two Shower Partons in a Jet in Vacuum 

When a hard parton i with any fixed large momentum k fragments in vacuum to a shower 
of semi-hard and soft partons, we denote the single- and two-parton distributions in terms 
of the momentum fractions X\ and X 2 by 

Pi( 1 ) = ^i(w) , (17) 

P 2 (l, 2) = (ssy/(x 1 ,x 2 ) , (18) 

where the RHS of Eq. (18) is to be identified with Eq. (7). For definiteness we consider 
j = u and j' = d. For % — g we get for r 2 (l, 2) a 2D distribution as shown in Fig. 4. Other 
distributions with i = u or d are very similar. 

Since r 2 (l,2) = 1 implies no correlation, we see from Fig. 4 that the two shower partons 
are uncorrelated when x\ ~ X 2 ~ 0, but are highly correlated at higher x t . Kinematical 
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constraint cine to momentum conservation forces r 2 (l,2) to vanish at x\ + x 2 — 1. It 
approaches zero faster, as x\ = X 2 —■> 0.5, than it does when X 2 is fixed and small, as 
X\ —> 1 — X2- The implication is that if these two partons are to recombine to form a pion 
at x = X\ + X2, there is more contribution from the asymmetric pairing of X\ and X2 than 
from the unnecessarily stringent requirement that x\ = X2 = x/2. 

3.2 Two Shower Partons in a Jet in HIC 

For each event in HIC with a high px jet, we can calculate the two-parton distribution as 
above. However, the hard parton momentum k cannot be fixed in an experiment. The distri¬ 
bution in the shower parton momenta q\ and g 2 , when averaged over all events, necessitates 
an integration over k and a sum over all i. The singlc-parton distribution 

Mi ) = s j (qi) ( 19 ) 

is given by Eq. (5), while the two-parton distribution 

p 2 (l,2) = (55)"'(g 1 ,« ! ) (20) 

is given by Eqs. (6) and (7). Some details in the calculation need to be explained. When i 
in Eq. (6) is g, then SjSj = GG, but when i — u, say, then S{S{ = KL, where G, K and 
L are SPDs initiated by gluon, valence quark and sea quark, respectively, in the notation of 
[10]. It turns out that the gross features of the results do not depend sensitively on these 
details. 

For Au+Au collisions at 200 GeV we have calculated r 2 (l, 2). The results for most central 
and most peripheral collisions are shown in Fig. 5. The main feature is that r 2 (l,2) is large 
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and increases with q\ and g 2 . That is markedly different from Fig. 4. The reason is that 
Pi(l)pi(2) depends quadratically on £/*(&), whereas p 2 (l,2) depends on it linearly. The 
probability £_/)(&) for a hard parton to be produced at k in HIC is small in central collisions, 
and is even smaller in peripheral collisions (although £ is larger) because f t (k) scales with 
(TVcou). Thus the probability of producing two partons in a jet is much higher than the 
product of the probabilities of producing the two partons separately. The latter does not 
require the two partons to be produced in two jets in the same event. For experiments at 
RHIC we may assume that events with two jets in the same hemisphere never occur, since 
they are so rare. Another way of interpreting r 2 (l, 2) is to write it as 

r 2 (l,2) = p c (2|l)/p 1 (2) , (21) 

where p c (2|l) is the conditional probability of finding a parton at g 2 given that a parton is 
at qi 

p c ( 2|1) = HIA . (22) 

PH 1 ) 

Thus p c (2|l) is like the distribution of a particle associated with a trigger at q±. The largeness 
of the ratio in Eq. (21) means that it is far more likely for a shower parton to be found at g 2 in 
association with another parton already found at q\ than to find a parton at g 2 independently. 

The two surfaces in Fig. 5 have essentially the same shape. The peripheral case is higher 
because the corresponding £_/)(&) is smaller. In Fig. 4 the hard parton momentum k is held 
fixed and the probability of having such a parton is 1, while in Fig. 5 k is integrated up to 
40 GeV/c so that at higher qi and g 2 there is no kinematical cut-off, albeit the probability 
fi(k) of producing a hard parton at higher k is greatly suppressed, which only makes r 2 (l, 2) 
higher. 
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The study in this section may seem academic, since shower partons are not directly 
observable. However, those shower partons are the main constituents that contribute to 
the formation of pions that we shall treat in the following section. Their properties are 
therefore important to help us understand the features that we shall find about dihardron 
correlations. It should be emphasized that no free parameters have been used in this section, 
or in the section to follow, apart from C, T, £ that have been determined in Sec. 2 for the 
single-particle distributions. 


4 Correlation of Pions in Jets 


For the production of two pions it is necessary to consider the recombination of four partons. 
The corresponding equation for the two-pion distribution is [9] 


dK 


7T17T2 


pip 2 dpidp 2 p\pl 


' 4 dq\ 

n — (? 2 , 93, q^Rmiqi, 93, Pi)R^ 3 {q 2 , P 2 ). 

«=i ?< J 


(23) 


The partons with momenta q\ and q 3 that form can be either thermal or shower, as in 
Eq. (3). Similarly, the same is true with q 2 and q 4 . Thus F 4 has the form 


F 4 = (TT + ST + SS) 13 (TT + ST + SS ) 


24 


(24) 


Among the various possible combinations of terms, it is clear that ( TT)i 3 (TT) 24 has nothing 
to do with hard scattering and contributes to the background. For that reason the terms 
involving (T T) are not included in [9], where the background-subtracted distributions of the 
associated particles are considered. Here we study C 2 (l,2), which has its own definition of 
subtraction, so we start out being general. Despite the deceptive appearance of Eq. (24), 
F 4 is not factorizable. If it were, its substitution into Eq. (23) would render the factorizable 
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form p-2^1,2) = pi(l)p2(2), implying no correlation. The nonfactorizability of F 4 is rooted 
in the correlations among the various shower partons. The correlation between two shower 
partons in (55) i 3 or ( 55)24 has already been considered in connection with Eq. (20), which 
is shown explicitly in Eqs. (6) and (7). However, those are not the correlations that render 
F 4 nonfactorizable. It is the correlation between the partons in (1,3) with the partons in 
(2,4) that gives rise to nonvanishing d7 2 (l, 2). More specifically, because all shower partons 
in Eq. (24) originate from the same hard parton at k , there is only one ffi(k ) factor for 
all 55 • ■ ■ terms. Thus, for example, in (5T) 13 (5T) 2 4 the (55)(gi,g 2 ) term is exactly as 
expressed in Eq. (6). In (5T) 13 (55) 2 4 we can append R W2 iq 2 , q 4 ,p 2 ) and make use of Eq. 
(8) to convert (SS) 2 aR-k 2 to the FF .D 71 " 2 , so that the (5T)(55) contributions to p2(l,2) is 
explicitly 


dN stss 

^ 7T17T2 


£ 


pip 2 dp\dp 2 P 1 P 2 


E [ dkkfiik) I dq 1 T(p 1 -q 1 ) 

A At Alt J J 


•E- 

i 2 


SI ( $ 


k ) k — qi 


-Dl) 


P2 


k - gi 


+ si 


qi 


k — P 2 ) k 


-DFf ( t- 
k 


(25) 


which is not factorizable. For definiteness we shall hereafter consider Ti\ and 712 to be both 
7r + . Then the sum in j is over u and d. If i = g, then i\i" — g\ if i — V (u or d), then 
= (V, 5) or (5, V) where S denotes sea; if i — 5, then i',i" = S also. 

For completeness, let us write out the two other combinations of terms. For (5T) 33 (5T)24 
we have 


d,N: 


stst 


pip 2 dp 1 dp 2 


P 1 P 2 w 


J 2 J dkkfi{k) J dqidq 2 T (pi - qf)T(jp 2 



92) 


(26) 
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where j and j' can be either u or d. For (SS) 13 ( 55)24 we have 


dN ssss 

** 7T17T2 




£ fdkkMk)l-D? 


Pi 


-D 


7T2 


P2 


(27) 


PiP^dpiipi PiP 2 jjrJ„J I 7 \kj k-p 1 ' \k-pij J 

In all three equations above i, i' and i" follow the same rule stated at the end of the previous 
paragraph. For simplicity let us use the symbolic form for these three equations as follows: 
[(ST)T>\, [(«ST)(«ST)] and [DT>], respectively. The parentheses imply the contents are to 
recombine; the square brackets signify that the shower partons inside all arise from one 
hard-parton jet. 

It is appropriate at this point to compare our description of two-particle distributions 
with the correlation study of Fries et al (FBM) [17]. Since both approaches are based on the 
recombination model [6, 18, 19], there are many similarities, particularly in the separation 
into various components. Apart from the difference in focus where we consider correlation in 
Pt, while FBM studies correlation in 0, the main difference in physics is that we break the 
hard parton into shower partons, whereas FBM does not. If our TS recombination is related 
qualitatively to their soft-hard (SH) recombination, then the classification of components is 
essentially the same in the two approaches. Quantitatively, TS can be significantly different 
from SH. Moreover, the correlation at the parton level is for us among the shower partons, 
while for FBM among the soft partons. 

In /? 2 (1, 2) the (TT){TT) term is factorizable, and is cancelled by the thermal terms in 
pi(l)pi(2), where pi(l) is given by Eq. (10). That corresponds to background subtraction. 
Other double TT terms like ( TT)(ST) also get cancelled. The surviving terms in C*2(l,2) 
all involve shower partons in each pion. We group them in four clusters of terms 


C 2 (l, 2) = Cf st ( 1,2) + Cf ss ( 1,2) + C s 2 sst { 1,2) + C| MS (1,2) , 


(28) 
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where 


C|“(l, 2) = [(ST)(ST)] 12 - (ST), » (ST), (29) 

C?"( 1,2) = [(Sr)(C )] 12 - (ST), » (C ) 2 (30) 

CT‘(1,2) = [D(ST)] 12 - (®)i » ( STh (31) 

Cr(l,2)=[»B] 11 -(B) 1 *(B) 2 • (32) 


The products involving * signify terms from p 1 (l)p 1 (2), so their shower partons have their 
own hard partons. More precisely, (ST )i * (ST) 2 , for example, has the explicit expression 

(ST), » (STh = (AY. j dkkm jdq,S’ T(p, - ?,) 

Jr £ / dk'k'Uk') I dq 2 S{ T(p 2 - q 2 ) j . (33) 

Clearly, this is very different from [( ST)(ST )] l2 , written out in full in Eq. (26). Thus 
Cf st ( 1,2) does not vanish, in general, and likewise for the other terms in Eq. (28). 

Regarding Eqs. (27) and (32) that involve FFs only, it is relevant to compare them to the 
medium modihed dihadron FF studied in [20]. In our formulation of hadron production in 
HIC the medium effect is taken into account by considering (TT) and (TS) recombination, 
leaving the (SS) recombination to be unaffected by the medium except for the hard parton 
distribution fi(k). Thus in Eqs. (9) and (27) the FF in vacuum is used. The contribution of 
C'| sss (l, 2) to the observed correlation cannot be measured separately, except in pp collisions 
or in very peripheral AA collisions where the medium effect is negligible. In the fragmentation 
model modification of the FF is the only way to take into account of the medium effect, 
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whereas in the recombination model fragmentation in vacuum is unimportant compared 
to thermal-shower recombination. For these reasons it is inappropriate to compared Eq. 
(27) for central nuclear collisions with the medium modified dihadron fragmentation process 
considered in [20] without including Eqs. (25) and (26) also. 

We have calculated all the quantities in Eqs. (28) - (32). The results cannot easily be 
displayed because each term decreases rapidly with p\ and p 2 and the various clusters of 
terms, C%, {a = stst ,..can be positive and negative, so that neither linear nor log plots 
are appropriate. To give a sense of the complication, we show in Fig. 6(a) the result for 7r + 7r + 
correlation in central Au+Au collisions for the restricted case of pi = p 2 = p in the range 
1.5-10 GeV/c. At p ~ 2 GeV/c, all Cf become negative, resulting in a dip in C 2 . There is 
no such fine structure for peripheral collisions (80-92%), so a semi-log plot can exhibit the 
behaviors of C% in the entire range up to p — 10 GeV/c, as shown in Fig. 6(b). Note that for 
0-10% centrality the terms involving st component is dominant, while for 80-92% centrality 
the a = ssss component is dominant. These features are in accord with the properties 
of single-particle distributions shown in Fig. 3, where thermal-shower and shower-shower 
components, respectively, are dominant in the two cases. 

When pi and p 2 are different, the two surfaces of C 2 (l,2 ) for the two centralities are 
shown in Fig. 7 in the 2D ranges where they are positive. The region where C*2(l,2) is 
negative cannot be shown in the semi-log plot. For what can be shown they behave rather 
similarly. They both are suppressed by about 10 orders of magnitude going from 1 to 10 
GeV/c. A large part of the rapid decrease is due to the power-law suppression that is present 
even in the single-particle distributions over the same Pt range, as can be seen in Fig. 1, 
though not more than 8 orders of magnitude. Thus the correlation function C2(l,2) is not 
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the best quantity to exhibit the properties of correlation. 

To de-emphasize the power-law suppression of the hard parton distribution at high k 
that is contained in pi(l) and pi(2), we consider the normalized correlation K 2 ( 1,2) defined 
in Eq. (14). When we set Pi = p 2 = P, the dependence of K 2 (p) on p is shown in Fig. 8 for 
(a) 0-10% and (b) 80-92% centrality. In (a) a dip occurs at low p, the minimum being at 
p = 3 GeV/c with K 2 < 0; in (b) there is no dip. The rapid rise at larger p is introduced 
into K 2 ( 1,2) by the division of C 2 (l,2) by pi(l)pi(2), which is strongly damped at high pi 
and p 2 . It is the result of an over-correction. The fact that A^SO — 92%) S> AT 2 (0 — 10%) 
at all p in Fig. 8 is due to the property of £ fi(k ) being much smaller in peripheral collisions 
compared to central collisions. Since the numerator of K 2 ( 1,2) depends linearly on £ fi{k ), 
but the denominator depends quadratically on it, the peripheral K 2 ( 1,2) is therefore much 
larger. 

What we see in Fig. 8 (a) and (b) about the general features of AT 2 (1,2) for the two 
centralities is very similar to what is shown in Fig. 5 for r 2 (l, 2) for shower partons, recalling 
that A' 2 (l, 2) = r 2 (l, 2) — 1, as stated in Eq. (14). Although r 2 (l, 2) cannot become negative, 
it becomes large in Fig. 5 at high q\ and q 2 for both centralities. Of course, hadron correlation 
is far more complicated than parton correlation. But the common feature that they both 
become large at large parton and hadron momenta is due to the commonality in the definition 
of r 2 (l,2) and A' 2 (l,2), which have pi(l)pi(2) in the denominators. 

It is now clear that we need something that is between C 2 ( 1, 2) and A' 2 (l, 2) so as to have 
linear dependence on £ f^k ) in both the numerator and the denominator. That quantity is 
the ratio G 2 ( 1, 2), defined in Eq. (16). Figure 8 (c) and (d) shows the dependence of G 2 (l, 2) 
on Pi = p 2 = p for (c) 0-10% and (d) 80-92% centrality. Now, G 2 (p) is small at large p in 
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both cases, and the properties of both can be revealed in linear plots. For 0-10% centrality 
there is still a dip at p — 2 GeV/c with £ 2 (p) < 0 in its vicinity, as it should, but for 80-92% 
centrality G 2 (p) decreases monotonically with increasing p. The normalizations of G 2 (p) for 
the two cases are now nearly the same. That can be seen more easily in semi-log plots. In 
Fig. 9 we show all three correlation function, C 2) K 2 and G 2 , for the two centralities. The 
gap in Fig. 9(a) is due to the functions being negative in that range. Evidently, we have 
£2(0-10%) » £2(80-92%), W 2 (0-10%) < Jl 2 (80—92%), but £ 2 (0-10%) « £ 2 (80-92%) 
in regions where they can be compared. Thus the function £ 2 (1,2) is best for the study of 
the nature of correlation. 

In Fig. 10(a) we show £ 2 (1, 2) in 2D distributions for both centralities in the same linear 
space. We now can see clearly their differences, knowing that asymptotically they behave 
similarly without structure. In the ranges of pi and p 2 plotted £ 2 (1,2) shows negative 
correlation for 0-10%, but positive correlations for 80-92% centrality. That difference can be 
succinctly exhibited in the ratio 


Rcp( 1 . 2 ) = 


G<o-io%>(1,2) 


^,( 80 - 92 %) 

^2 


( 1 , 2 ) 


(34) 


which is shown in Fig. 10(b). Since peripheral Au+Au collisions at 80-92% centrality is very 
close to pp collisions, we may regard Fig. 10(b) as the clearest revelation of the difference in 
dihadron correlation between central Au+Au and pp collisions. 

The origin of the dip in Fig. 10(b) can be found in Fig. 6 where £ 2 is negative at 
Pi = P 2 = P = 2 GeV/c for 0-10% centrality, but not for 80-92% centrality. Let us now 
describe the reasons why C 2 are negative in Fig. 6(a), but not in Fig. 6(b). Many factors 
are involved with competing strengths that can be compared only by concrete calculation, 
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but qualitative discussion can provide some insight. Consider Eq. (29) whose two terms are 
shown explicitly in Eqs. (26) and (33). The former has one factor of £fi(k) and two SPDs 
sharing the same k, while the latter has two independent factors of £fi(k), each having only 
one SPD. In the peripheral case £fi(k) is very small because ,f){k) is scaled down by (1 V co ii), so 
{ST )i * {ST) 2 is always small compared to [( ST){ST)\ V2 at all qy and q 2 . Hence, C' 2 stsi (l, 2) 
is always positive; the same argument applies to all other (7^(1, 2). In the case of central 
collisions £ fi{k ) becomes large enough to allow other factors to also play some role. When 
Pi and p 2 are on the lower end of the intermediate range, e.g. < 4 GeV/c, qi and q 2 must 
be even less, but not too small because T{p x — qk) and T{p 2 — q 2 ) would be exponentially 
damped at large pi — q % . In other words, if the momentum fraction ay of the shower partons 
are small, there are not many thermal partons to make up p\ and p 2 \ on the other hand, if ay 
are large, the SPDs are suppressed. The hard parton k cannot be large to relieve the latter 
constraint without suffering the suppression of fi{k). This situation of multiple constraints 
is easier for {ST)\ * {ST) 2 to meet since there is only one shower parton per jet in Eq. 
(33), but is much more difficult for [{ST){ST)\i 2 where the jet momentum must be shared 
by two shower partons. Thus [{ST){ST)\ 12 is forced to be smaller than {ST)i * {ST) 2 in 
central collisions where kfi{k) is not infinitesimal and thermal partons are important in the 
intermediate Pt region. 

In the case of (7f sss (l,2) there is no T, but the demand for additional shower partons 
plays the same role. It is clear from Eq. (32) that double fragmentation causes \DT>]i 2 to 
be more suppressed than two single fragmentation {V )i * {T>) 2 . Thus (7f sss ( 1,2) becomes 
negative so long as p\ and p 2 are not high enough to force the large k to dampen double £,fi{k) 
in {T >)i * {V) 2 more severely than the single £fi{k) in [DT>]i 2 . The other C% terms behave 
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similarly for the same reason, resulting in the sum 6 * 2 ( 1 , 2 ) to be negative when p\ and p -2 are 
both less than around 3 GeV/c. The negativity of 6 * 2 ( 1 , 2) should not be surprising, since 
the ratio r 2 (l, 2) for the shower partons in a jet of fixed k is less than 1, as shown in Fig. 4. 
The corresponding 6 * 2 ( 1 , 2) is negative for all X\ and X 2 ■ In the hadronic problem there are 
in some terms even more shower partons involved, all competing for the limited momentum 
that the hard parton provides. Fig. 10 exhibits the intricacies of the hadronic correlation 
even when the basic correlations among the shower partons are kinematical in nature. 


5 Comparison with Experimental Data 

In the foregoing we have studied the correlation of two 7 r + without using any adjustable pa¬ 
rameters. What we have obtained unfortunately cannot be compared directly with current 
data, because most experiments either do not yet have particle identification, or put empha¬ 
sis on jets with trigger particles and study the dependencies on A<p and A// of associated 
particles. There are, however, two areas in which our theoretical calculations come close to 
the experimental data and are of interest for us to examine them in more detail here. 

The first area concerns the experiment in [5], in which data are given for the px distri¬ 
bution of charged particles associated with trigger in the range 4 < p T < 6 GeV/c. We can 
come close to what is measured by calculating 


dN‘ 


tr—bg 


f 6 , dN*+Jl . r 6 , dlW 

= / dp ip i - T / / dpipi 

J 4 PlV2dPldV2 2r 


(35) 


p 2 dp 2 h ^p\P2dp\dp2 h 'pidpi 
where p\ is regarded as the momentum of the trigger particle, and P 2 that of the associated 
particle. The superscript, tr — bg, denotes the sample of triggered events with background 
subtracted. For the integrand in the numerator we use Eq. (23), but with F 4 not given 
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by Eq. (24). Background subtraction corresponds in our case to leaving out in Eq. (24) the 
terms that are faetorizable. We are uncertain whether that correspondence is exact, since the 
experimental procedure for the determination of the background is not unique and involves 
steps that seem hard to duplicate in theoretical calculation. Onr procedure of retaining only 
the non-factorizable component is the most sensible way to obtain the associated particle 
distribution, so we use 

ptr-hg = {sr + SS) 13 (ST + SS) 24 , (36) 

which is not faetorizable because of the correlation that exists between the shower parton at 
qi and the shower parton at q 2 , as we have studied in Sec. 3. 

We note that Eq. (35) is a ratio of the integrals over the trigger momentum in contrast 
to the integral of the ratio that is calculated in [9]. Eq. (35) corresponds to the way that the 
data are presented in [5]. 

There is one piece of detail about the calculation using Eqs. (23) and (36) that we now 
discuss. Among the SPDs determined in [10] the gluon initiated quark distribution has a 
parametrization that leads to divergence, G(x) oc x~ a , as the momentum fraction x —> 0. ft is 
unphysical in that limit, although the parametrization renders a good fit of the fragmentation 
function away from that limit. In application to the calculation of dN n /pdp, as in Eq. (10), 
the integral over q\ is convergent, and at small p the ST component is dominated by the 
TT component. But when the background is subtracted with the (TT)(TT) component 
removed from dN 7Tn /pip 2 dpidp 2 , the effect of large G(x) at very small x is exposed. Since 
that effect is unphysical and since we do not claim reliability of our formalism at very small 
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Pt, we introduce an ad hoc factor 


s(q) = 1 - e~ q/qo (37) 

to suppress the distribution at small q for all gluon initiated shower partons. With q 0 set at 
1 GeV/c this procedure does not affect the results of our calculation for pt > 1.5 GeV/c, 
but does introduce a round-off at lower p T in the associated particle distribution. Figures in 
Sec. 4 all have that factor s(q) included for gluon initiated showers. 

Substituting Eq. (36) into Eq. (23) and then into Eq. (35), we obtain the results shown 
in Fig. 11 for (a) 0-10% and (b) 80-92% centralities. The calculated curves are for n + trigger 
and 7r + associated particle, while the data points from [5] are for all charged particles with 
0-5% centrality in (a) and for pp collision in (b). The shapes of the Pt dependence agree 
very well; the normalizations differ by a factor that can be attributed to the data being for 
all charged particles. Thus our formalism describes the correlation between particles in jets 
quite satisfactorily when compared to data. 

Another area of possible comparison between theory and experiment is in correlation at 
low pt- Data on a quantity denoted by A p / v /p in Ref. [2] for pp collision show a bulge 
plotted in transverse rapidity, defined by 

y t = ln[(m T + p T )/m n \. (38) 

Since the quantity Ap / v /p is rather similar to our G^l, 2), it is of interest to display the latter 
in terms of y t , which has the virtue of expanding the low pt scale. Since our calculations are 
unreliable for pt < 1 GeV/c, we consider only the region y t > 2.5. Instead of Eq. (13), we 
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now define 


Pi(l) 


dN wl 

d'l-Ji 


p 2 ( 1,2) 


(l 1^-Jl i TV 2 

dyidy 2 ’ 


(39) 


where dyi = dpi/rriT■ C 2 ( 1,2) and G 2 ( 1 , 2 ) remain to be defined by Eqs. (12) and (16), 
respectively, but in terms of y\ and y 2 they are changed in accordance to Eq. (39). 

In Fig. 12 we show G 2 (l,2) as functions of y x and y 2 for (a) 80-92% and (b) 0-10% 
centralities. They are reproductions of the two surfaces in Fig. 10(a), now in terms of y,. 
The upper ones in both are for peripheral collisions, while the lower ones are for central 
collisions. The increase at small y\ and y 2 in Fig. 12(a) appears to agree with the bulge in 
the pp collision data in [2] for yi,y 2 > 2.5. For central collision Fig. 12(b) has a peak and a 
dip, which are predictions that should be checked by the analysis of the corresponding data. 
For completeness, we show in Fig. 13 the R^p in the yi,y 2 variables, in terms of which the 
dip becomes the dominant feature, even more prominent than in Fig. 10 (b). 


6 Conclusion 

We have calculated the single-particle distribution in Au+Au collisions at all centralities, 
and considered various types of correlation functions in jets. For correlations between shower 
partons we have shown the drastic difference between the case when the hard parton mo¬ 
mentum k is fixed, as in e + e~ annihilation, and the case in heavy-ion collisions where k is 
not fixed. For hadron correlations we have considered three types of correlation fucntions, 
(72(1,2), K 2 (l,2), and G 2 (l,2), and showed that G 2 ( 1,2) is best in displaying the structure 
of the dihadron correlation. 


23 



We have found that the ratio of 6*2(1, 2 ) for central collisions to that for peripheral colli¬ 
sions has a big dip in the p\ ph p-2 ~ 2 GeV/c region. That is where the negative correlation 
between shower partons in a jet pulls down the hadronic p 2 (l, 2 ) in central An+Au collisions 
relative to pi(l)pi( 2 ). However, in peripheral collisions pi(l)pi( 2 ) is so small compared to 
P2(l,2) that 6*2(1, 2 ) remains positive for all p\ and P2. The calculated features of 6*2(1, 2 ) 
can be checked by experiment, since all the quantities involved are directly measurable. 

Although most of the results obtained on the correlation functions cannot be compared 
with any existing data, there is no basic hurdle to overcome for the data analysis to be 
carried out in ways that can check onr prediction. The distribution of associated particles 
that we have calculated for 7r + 7r + correlation compares satisfactorily with the data except for 
the normalization difference dne to the fact that the data are for all charged particles. Onr 
result on G2 for peripheral collisions has some resemblance to the data on pp collisions, which 
are, however, not presented in a way that can facilitate quantitative comparison. Clearly, 
closer coordination between theoretical and experimental efforts would be of great value to 
both. 

In this work the basic quantity is the correlation function 6*2(1, 2) that is precisely de¬ 
fined. By subtracting pi(l)pi( 2 ) from p 2 (l, 2 ), not only is the background subtracted, other 
factorizable terms are taken out too, creating the possibility that 6*2(1, 2) can become neg¬ 
ative. That possibility does not show up in the present experimental way of determining 
the associated particle distribution [ 3 , 5 ], but may be present in other ways of studying 
correlations [ 1 , 2 , 4 ], We suggest that future experimental work be directed toward deter¬ 
mining the correlation function 6 ^( 1 , 2 ) and the central-to-peripheral ratio Rcp{ 1 , 2 ). The 
absence of any hole in the experimental result on Rcp{ 1 , 2 ), as shown in Figs. 10 (b) and 
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13, would provide valuable information on the nature of parton correlation that we have not 


incorporated in our formulation of the problem. 


Acknowledgment 

We are grateful to Tom Trainor and Fuqiang Wang for valuable communication on their 
work and to C. B. Yang for helpful comments. This work was supported, in part, by the U. 
S. Department of Energy under Grant No. DE-FG03-96ER40972. 

References 

[1] T. A. Trainor (STAR Collaboration), talk given at the 20th Winter Workshop on Nu¬ 
clear Dynamics, Jamaica, 2004, hep-ph/0406116. 

[2] R. J. Porter and T. A. Trainor, Acta Physica Polonica B 36, 353 (2005). 

[3] D. Magestro, Proceedings of Hard Probes 04, Ericeira, Portugal (2004). 

[4] J. Adams et. ah, (STAR Collaboration), nucl-ex/0406035, 0408012, 0411003. 

[5] J. Adams et. ah, (STAR Collaboration), nucl-ex/0501016. 

[6] R. C. Hwa and C. B. Yang, Phys. Rev. C 70, 024905 (04). 

[7] R. C. Hwa and C. B. Yang, Phys. Rev. Lett. 93, 082302 (04); Phys. Rev. C 70, 037901 
(04). 

[8] R„ C. Hwa, C. B. Yang and R„ J. Fries, Phys. Rev. C 71, 024902 (2005). 


25 



[9] R, C. Hwa and C. B. Yang, Phys. Rev. C 70, 054902 (2004). 


[10] R„ C. Hwa and C. B. Yang, Phys. Rev. C 70, 024904 (2004). 

[11] R. C. Hwa and C. B. Yang, hep-ph/0312271vl. 

[12] B. A. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B597, 337 (2001). 

[13] R. C. Hwa and C. B. Yang, Phys. Rev. C 66, 025205 (2002). 

[14] D. K. Srivastava, C. Gale, and R. J. Fries, Phys. Rev. C 67, 034903 (2003). 

[15] S. S. Adler et al (PHENIX Collaboration), Phys. Rev. Lett. 91 , 072301 (2003). 

[16] E. W. Weinsstein, Correlation Coefficient, MathWorld - A Wolfram Web Resource 
http: / / mathworld.wolfram.com/CorrelationCoefhcient.html. 

[17] R. J. Fries, S. A. Bass, and B. Miiller, Phys. Rev. Lett. 94, 122301 (2005). 

[18] R. C. Hwa and C. B. Yang, Phys. Rev. C 67, 034902 (2003). 

[19] R. J. Fries, B. Miiller, C. Nonaka, and S. A. Bass, Phys. Rev. Lett. 90, 202303 (2003) 
Phys. Rev. C 68, 044902 (2003). 

[20] A. Majumder, E. Wang and X. N. Wang, nucl-th/0412061. 


26 



Figure Captions 


Fig. 1. Single-particle distributions of produced n in Au+Au collisions at •>fs = 200 GeV for 
four centrality cuts. Data are from [15], and the solid lines are the results of calculation 
in the recombination model. 

Fig. 2. The three parameters C, T and £, determined by fitting the data at all centralities. 

Fig. 3. The single-particle distributions for (a) most central and (b) most peripheral colli¬ 
sions, showing the three components that contribute to the sums in solid lines. 

Fig. 4. (Color online) The ratio r 2 ( 1, 2) in terms of the shower parton momentum fractions 
x\ and X 2 for a hard gluon momentum fixed at a large values k. 

Fig. 5. (Color online) The ratio r 2 (l,2) in terms of the shower parton momenta q\ and q 2 
in heavy-ion collisions for two extreme centralities. 

Fig. 6. The correlation functions C 2 for 7r + 7r + in HIC, showing the contributions from the 
three combinations of S and T partons for (a) central and (b) peripheral collisions 
along the diagonal where p\ = p 2 . 

Fig. 7. (Color online) The correlation functions C 2 for two centralities when they are 
positive. 

Fig. 8. The normalized correlation functions K 2 along the diagonal where pi = p 2 for (a) 
central and (b) peripheral collisions. The partially normalized correlation functions G 2 
for (c) central and (d) peripheral collisions. 
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Fig. 9. A comparison of the three types of correlation functions C' 2 , K 2 and G 2 in semi-log 
plots, with the gap corresponding to where the functions are negative. 

Fig. 10. (Color online) (a) G 2 (l,2) plotted as functions of p\ and p 2 for two centralities; 
(b) the ratio R^p of central to peripheral G 2 . 

Fig. 11. Associated particle distributions of tt + i r + for (a) central (0-10%) collisions com¬ 
pared to the STAR data [5] for all charged particles at (0-5%) centrality and (b) 
peripheral (80-92%) collisions compared to the STAR data [5] for all charged particles 
in pp collisions. 

Fig. 12. (Color online) G 2 ( 1,2) for 7r + 7r + plotted in terms of transverse rapidities for (a) 
peripheral (80-92%) and (b) central (0-10%) collisions. 


Fig. 13. (Color online) The ratio R^p in transverse rapidities. 
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